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Abstract 



We provide, in the R-parity violating super symmetric standard model, a 
comprehensive analysis for the determination of sneutrino vacuum expecta- 
tion values from the one-loop effective scalar potential, and also for one-loop 
renormalized neutrino masses and mixing by calculating the effective neutrino 
mass matrix in the weak basis. Applying our results to theories with gauge 
mediated supersymmetry breaking, we show how atmospheric and solar neu- 
trino oscillations can be accommodated simultaneously in this framework. It 
is observed that the one-loop correction to sneutrino vacuum expectation val- 
ues leads to a significant effect on the determination of the neutrino masses 
and mixing. 

PACS number(s): 12.60.Jv, 11.30.Fs, 14.60.Pq 
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I. INTRODUCTION 



The minimal supersymmetric standard model (MSSM) may allow for explicit lepton num- 
ber and thus R-parity violation through which neutrinos get nonzero masses and mixing [|TJ . 
As it is an attractive possibility to explain neutrino experiment data, there have been many 
works investigating neutrino masses from R-parity violation . In recent years, further de- 
velopment has been made to analyze three neutrino masses and mixing || accommodating 
the recent observation of the atmospheric muon neutrino oscillation in Super-Kamiokande 
0] and other experimental data A very interesting feature of R-parity violation as 

the origin of neutrino masses and mixing is that this idea could be tested in future collider 
experiments despite the small R-parity violating couplings. Specifically, the mixing angles 
measured in the Super-Kamiokande Q and the CHOOZ experiments can be probed by 
detecting lepton number violating decay modes of the lightest supersymmetric particle [|5],[| . 

In this paper, we basically extend the former studies in two ways. First, we calculate 
one-loop improved tree- level neutrino masses. In order to do that, we determine sneutrino 
vacuum expectation values (VEVs) from the minimization of the one-loop effective scalar 
potential, which has usually been ignored in previous works Q As will be shown later, the 
one-loop correction to sneutrino VEVs should not be ignored as it can drastically change 
the tree-level results. When there exists a cancellation among tree-level quantities in the 
determination of sneutrino VEVs, the loop contributions to the neutrino masses can be 
even larger than the tree ones. Furthermore, this phenomenon occurs generically in a certain 
parameter region. Secondly, we calculate one-loop renormalized neutrino masses and mixing 
utilizing the "effective mixing matrix" method |IJ . As far as one concerns only the neutrino 



sector, this approach turns out to be very useful even though the results should be the same 



as in the usual on-shell renormalization scheme working in the tree-level mass basis flTT| . |T2 
Our procedure begins with constructing the 7x7 one- loop neutrino/neutralino mass matrix 
in the weak basis. Then, we obtain the seesaw-reduced 3x3 neutrino mass matrix from 
which physical neutrino masses and mixing are calculated. In this way, one can easily figure 
out the effects of each self-energy diagrams and simplify the calculation to obtain the main 
contributions. Our results are applied to the theories with gauge mediated supersymmetry 
breaking to discuss how large the one-loop results deviate from the tree-level ones. Finally, 
we show that the atmospheric and solar neutrino oscillations can be accommodated naturally 
in this framework. In particular, we present some typical examples for input parameters 
providing the vacuum or matter oscillation solution to the solar neutrino problem. 

This paper is organized as follows. In Section II, we determine sneutrino VEVs by 
computing one-loop effective scalar potential, and present analytic expressions for all the 
one-loop contributions to sneutrino VEVs in Appendix A. In doing so, we also derive the 
diagonalization matrices for fermions or sfermions which get mixed further due to R-parity 
violation. In Section III, these diagonalization matrices will be used to calculate one-loop 
self-energy diagrams in the weak basis. This enables us to determine the dominant contri- 
butions to physical neutrino mass eigenvalues and mixing angles which could be relevant 



1 see however, Ref. [12] where diagrammatic method is used. 
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for explaining the current neutrino data. The leading self-energy diagrams are calculated 
in Appendix B. In Section IV, we apply the above results to the R-parity violating theories 
with gauge mediated supersymmetry breaking and discuss the effects of one-loop corrections 
to sneutrino VEVs on the determination of the neutrino masses and mixing. Then, we illus- 
trate how the atmospheric and solar neutrino data can be explained simultaneously in this 
framework. We conclude in Section V. 



II. SNEUTRINO VEVS FROM ONE-LOOP EFFECTIVE SCALAR POTENTIAL 

We start by writing the MSSM superpotential Wq and W\ conserving and violating 
R-parity/lepton number, respectively: 

W = iiH^Hi + hlL&EZ + hfQiH^ + h^Q t H 2 Uf 

W 1 = ^UE 2 + Xij k LiLjEl + X'^LiQjDl . (1) 

An unpleasant feature of R-parity violation is that there are too many free parameters 
essentially uncontrollable unless a certain theoretical framework for Yukawa hierarchies is 
invoked. For later use, we assume that Aj = Aj33 and A- = A^ 33 dominate over the other 
couplings for the first two generations, guided by the quark and lepton Yukawa hierarchy. 
The bilinear model where no trilinear R-parity violating couplings are present at the fun- 
damental scale satisfies this assumption, as rotating away the bilinear terms generates the 
hierarchical trilinear couplings A, A'. 

One-loop effective scalar potential to determine the VEVs of the Higgs and sneutrino 
fields, (H® 2 ) and (L°) is given by 

V = \m\ iHi UE\ + BiLiH 2 + h.c] + m 2 Li \U\ 2 + V D + V 1 + ---, (2) 

where V D is the SU(2) x U(l) D-term, and V x = g^Str.M 4 (in ^ - §) is the one-loop 
correction. In our notation, rn 2 L . Hi contains not only the soft mass term but also the super- 
symmetric mass term /i/Zj from the superpotential (p]). Here we assume the R-parity violating 
parameters are small enough, /ij -C /i and rn1. Hi m|., etc., to produce light neutrinos 
of m u < leV, which would be a consequence of a horizontal symmetry Furthermore, 



we take the universality among soft supersymmetry breaking terms at the scale of super- 
symmetry breaking where the bilinear term /XjLji7 2 can be rotated away together with the 
corresponding soft terms. One of a natural framework for the universality is gauge mediated 
supersymmetry breaking which will be discussed in Section IV. Under these assumptions the 
intergenerational mixing terms like /ij/iiLjL] can be neglected. The minimization condition 
for the sneutrino fields, dV/dU} = 0, leads to 

/r0N _ B^H") +ml Hi (H°) + 
where E^ 1 are given by 



3 



(i) (2) gjj f , 



where | means to put L° to zero. This expression is valid when R-parity violating parameters 
are small enough, /ij <C /x, etc.. In the case of small R-parity violation, the Higgs VEVs can 
be obtained from the usual minimization condition neglecting R-parity violating couplings. 
The complete one-loop calculation for the Higgs VEVs can be found in Ref. Jl4| . 



R-parity violation mixes neutrinos, charged leptons, sleptons and sneutrinos with neu- 
tralinos, charginos, charged and neutral Higgses, respectively, and gives rise to more com- 
plicated and bigger mass matrices of fermions and sfermions. In the case of small R-parity 
violation, it is convenient to use first approximate seesaw diagonalization to rotate away the 
R-parity odd elements. Then one uses the further diagonalization process to obtain mass 
eigenvalues which contain the contributions from R-parity violating parameters up to their 
second order. It is now a simple manner to take derivatives of mass eigenvalues with respect 
to L° and find out the quantities in Eq. (|). 

Let us follow the above prescription first to get charged slepton/Higgs contribution to 
sneutrino VEVs. Since sleptons and charged Higgs bosons get mixed, one has to diagonalize 
the 8x8 mass matrix for three slepton pairs and one Higgs pair which takes the form: 

To calculate X^ 1 ' 2 ) we keep the R-parity violating contributions to the above mass matrix up 
to quadratic terms in R-parity violating parameters. In this paper, we assume the absence 
of intergenerational mixing in slepton soft masses. However, there could arise in general R- 
parity violating intergenerational mixing masses which can be neglected as they give quartic 
contribution to mass eigenvalues. We further assume CP conservation and thus take the 
VEVs, Ui = Vi t 2 = (Hi 2 ) to be real. Then, the submatrices of Eq. in the weak 

basis (e iL , e iR , H{, H%) are given by 

u2 = f M 2 Ll - (gl + \g 2 z )u 2 + glut - 2h\v 1 (\ i u i ) M 2 m + S^(A jUj ) 

11 \ Ml, + 5 i3 (A jUj ) Ml + (gl - gl)v? + (h\ Ul f - 2h%Vi(\ jUj ) 

M 2 - =( m ^ + ^" ~ h f) VlU * ~ Bi + & UiV2 ) (6) 
~ lH V hfv 2 Hi + A\ui hfvxfii - hf/xui J 

.. 2 (M 2 H --±(2g 2 w -g 2 )u 2 + h 2 T u 2 -B + gfao* \ 

HH ~\ ~B + glviv 2 M 2 Ht + \(2gl - g 2 )u 2 ) ' 

where Ml, etc are the conventional R-parity conserving charged slepton or Higgs masses 
and summing over j is understood for the quantities inside the parentheses in the first 
matrix, and we will take the lepton Yukawa couplings, h\ 2 = 0. Here we define the notation 
g 2 w = M^/v 2 and g\ = M\jv 2 where v 2 = v\ + v \ + u 2 and u 2 = J^i^ 2 - Following 
the previous prescription, one can derive the charged slepton and Higgs contributions to 
X^ 1 ' 2 ) as presented in Appendix A. Here we have to mention a few technical points. For 
diagonalizing away the R-parity odd matrix elements, we use the seesaw formula. That is, 
the matrix M. of the form; 
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KA f M l A 



can be diagonalized by the "approximate" unitary matrix U; 



U 



-±eet 



-e 



i-i©t©/ 



(7) 



(8) 



if M 12 and A are submatrices satisfying A <C M 12 . Being valid up to quadratic order in 
A, the seesaw rotation matrix is determined by the relation, A = MiG — @M 2 and the 
block diagonal matrices become 



Mi 
M 2 



M x + - (0A f + Ae f ) 
M 2 - | (0t A + A^) 



(9) 



The essential step in the diagonalization is to decouple the Goldstone mode, G~ . This can 
be done by first going to the basis of vanishing sneutrino VEVs with the rotation matrix Bg 
and then applying the usual rotation Up with the angle (3. Omitting the quadratic terms in 
R-parity violating parameters for simplicity, the combined rotation QcUp is given by 



/ 1 a t 

10 

-a, 1 

V 1/ 



1 



0\ 



-Cp Sp 

V0 s p cpj 



(10) 



where a, = Ui/v\ are the values determined from tree level scalar potential. Note that in the 

full expression, the angle (3 is determined by tan/3 = t> 2 /t4 where v[ = \jv\ + u 2 . Applying 
the above rotation to the charged lepton/Higgs mass matrix (^,^), one can explicitly see that 
there is a massless state up to the order of our expansion. Further diagonalization can be 
done by another seesaw rotation 0,5 which get rid of R-parity odd off-diagonal components, 
and the usual diagonalization U ei of the 2x2 slepton mass matrix [see Appendix A]. Again 
in the leading approximation, ®sU ei is given by 



QsU P 



(I 

1 

1 

\5\ 6i 





1 



s e 


V o 



Set 






0\ 



1 

1/ 



(11) 



We obtain the following expressions for 8\ 



2-' 



M 



EiH 



-aiSp(m 



Ml) - \ 

cp 



M% CH 



a,- 



fl Cp 



{m\ x -m\_){m\ a -m\S) 



8\ 



M% H (Ml-ml-)-Ml iH Ml t 



(ml u - m\_ 



){ml-m 2 H .) 



(12) 



In sum, the rotation matrix which brings the weak eigenstates (en, em, H 1 ,H 2 ) to the mass 
eigenstates (en,ei2,H~) decoupling the Goldstone mode G~ is 
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u c = 



-8\ + aiS[3\ 



V 



-Oi + 5\Sp)C ex - 5 t 2 S( 5 S ei 



-en + 5\s p )s ei + 5\s p c ei 
5\cps ei + 5 l 2 cpc ei 



-5\ 

C /3 



(13) 



This matrix will be used to compute one-loop self-energy diagrams for neutrino/neutralino 
masses. 

Following the similar steps, we also obtain the CP-even and CP-odd sneutrino/Higgs 
contributions to sneutrino VEVs. The mass matrix in the weak basis (z/f, iff*, iff) of the 
real components of slepton and Higgs fields is 



( Ml + &u 2 + g 2 z uj m 2 LiHl +g 2 z v lUi B t 



M 



u R H 



m 2 LiHl + g 2 z v lUi M 2 H n 



+ 4u 2 



B 



\ 



g 2 UiV 2 



B 



g 2 z vtv 2 



Ml* 



■ g 2 z UiV2 

in 2 



(14) 



where M 2 , etc are the usual R-parity conserving sneutrino/neutral Higgs boson masses. The 
diagonalization matrix in the leading expansion in terms of R-parity violating parameters is 



^ 1 Cl s a C 2 C a Cl c a C2 S a^ 

Cl S a Ca 

\ C2 C a S ct J 



(15) 



where a is the usual neutral Higgs scalar diagonalization angle [see Appendix A] and the 
seesaw rotation angle C\ 2 is given by 



c 2 = 



-<H(™% + M IWa4 - B M m % + M l C 2/3) 



m % ~ ™l. (ml + m 2 H )+ m 2 h m 2 H 
ai (m\ - M 2 c 2p )m 2 A s w /2 + Bi(m% - M 2 c w ) 



ml.(m\ 



m 



11) 



■mfm 



(16) 



h" l H 



As for the CP-odd sneutrino/Higgs fields, the mass matrix in the weak basis (z/ 7 , iff, H 2 ) 
of the imaginary components of slepton and Higgs fields is 



(Ml + iu 2 



M D i H = 



m 



V 



Li if! 

-Bi 



m L iHl 

M 2 Hf + Slu 2 
-B 



-B t 
-B 

M 2 r ~ k 



-w- 



\ 



(17) 



where M?i, etc are the usual R-parity conserving sneutrino/neutral Higgs masses. As in 

i 

the charged slepton/Higgs case, we first make the seesaw rotation, composed of = Uijv\ 
to decouple the Goldstone mode, G°. The diagonalization matrix which brings the weak 
eigenstates to the mass eigenstates (z> 7 , A ) becomes 



1 -Q + aiSp' 

U 1 ' = \ -Oi + Qsp s p 

<3 c /3 Cp 



(18) 
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where the seesaw rotation angle Q is given by 



C3 = 7 2 2~T • v ") 

In Appendix A, we collect all of the information needed to calculate one-loop corrected 
sneutrino VEVs and present complete analytic results. For the numerical calculation, we 
take the renormalization scale, Q 2 = m^m^, around which the one-loop contribution V\ is 
minimized. The mixing matrices given in Eqs. ([13]), ( fL5|) and (|T8|) play a role of inducing 
the effective A or A' couplings in the mass basis and will be used to calculate the one-loop 
neutrino masses. We also need the neutrino/neutralino and charged lepton/chargino mixing 
matrices |15,E| which we summarize in Appendix A for completeness. 



III. ONE-LOOP CORRECTED NEUTRINO MASS MATRIX 

The neutrino/neutralino mass matrix at one- loop level is given by 

e (p 2 ) = M{Q) + U(p 2 ) - i (M(Q)S(p 2 ) + £(p 2 )M(Q)) (20) 



where Q is the renormalization scale, and M(Q) is the the DR renormalized 7x7 tree-level 
mass matrix, II and £ are the contributions from one-loop self-energy diagrams. To compute 
the one-loop corrected mass matrix ( p0|) , we work in the weak basis, that is, we evaluate 
ITs and E's with weak eigenstate neutrino/neutralino on external legs and mass eigenstates 
running inside loops. For this, we express the propagators in the loops in terms of the mass 
eigenstates as calculated in Appendix B. A complete one-loop analysis has been performed 



in Refs. [11,12] where the mass matrix (2C) is calculated in the tree- level mass basis and mass 
eigenvalues and mixing matrix are obtained by re-diagonalizing the one-loop mass matrix. 
In our approach, it is useful to make a few observations which facilitate out calculation. Let 
us recast the mass matrix (|20T) as 



«"V> = G& ^) W) + (nf n> 



1 ( M D Y? D + Y> D Ml Y, n M D + M D T, N + Y, D M N . , ., 



N 



where neutrino-neutralino mixing elements of Mo or I1£),£,d are much smaller than M/v 
or ilAr,S7v by factor of 0(ai). If we are interested only in neutrino mass eigenvalues and 
mixing matrix, it is convenient to obtain first the seesaw reduced 3x3 neutrino mass matrix. 



Applying seesaw formula to Eq. ( 21[) and keeping one-loop order terms, we find the effective 



3x3 neutrino mass matrix as follows; 



MV) = -M D M^Ml(Q) + M D M^Ii^p 2 )M^ l Ml 

+ (M D M^M££ n (p 2 ) + E^MnM^Ml) (22) 
+ (M D M^ 1 Il£(p 2 ) + H^M^Ml) + n n (p 2 ) 
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As neutrino masses are negligible, we can now take p 2 = without ambiguity to calculate 
physical neutrino masses and mixing. Our procedure is in fact identical to that developed 
in Ref. |10[ where the mixing renormalization is performed by utilizing the Q-independent 



"effective mixing matrix". The effective mixing matrix U e ^ is calculated from the diag- 
onalization matrix U{p 2 ) of the mass matrix M pole (p 2 ) in the weak basis ( p0|) by taking 
Uij = Uij(m 2 ). According to the results of [|Kj, this method has several nice features. 
First, one can avoid superficial singularity reflecting the arbitrariness in the diagonaliza- 
tion of the mass matrix with degenerate massless neutrinos at tree level. Second, as far 
as neutrino mixing is concerned, the mixing matrix U u = U(p 2 —>■ 0) determined from the 
diagonalization of the effective 3x3 matrix ( p^) can be considered as the full unitary mixing 
matrix, up to a negligible correction of 0(a 2 ). Furthermore, the scale dependence which 
appears only in S n turns out to be very small as we will discuss below. 

Note that the first term in Eq. (|22|) is the tree-level mass which makes only one neutrino 
massive. The other terms show one-loop contributions from various self-energy diagrams. 
It is interesting to see that Sz^tv do not appear in M v {p 2 ) and thus one does not have to 
calculate them. We can make further simplifications. First, in the basis where /ij terms 
are rotated away, not only the tree-level mass matrix but also the loop mass matrix in 
Eq. (^) take a factorized form oc a^aj ||, That is, the contribution from II at is aligned 
with the tree-level mass matrix and thus gives overall correction to the dominant tree-level 
mass in the direction ~ a^Ui. One-loop correction to neutralino masses has been analyzed 



in Ref. |16|], according to which the II jy contribution typically gives 6% correction to the 
tree-level values. As we are interested in the leading correction to the masses of the two 
light neutrinos, we will ignore this correction by working in the basis with /Xj = 0. Secondly, 
once no flavor changing sfermion masses are assumed, S n has only diagonal elements which 
produces also aligned correction of the order a 2 jAm < 1%, and thus we can also ignore this 
contribution. Then we are left with the last two one- loop corrections involving U n and ILj 
in Eq. (^) which will give rise to a sizable contribution to neutrino masses and mixing. 
Therefore, the resulting neutrino mass matrix is 

M 2 M 2 
M% = + -j^&Vj + Vi&cp + U ViUj where 

F N = -{M X M 2 /M^ + M%s 2f >/n), Mj = M lC 2 w + M 2 s 2 w (23) 
-M 2 s 2 v M lC 2 w sp -c/3 

Note that £j = Oj — \i becomes Oj in the basis with /ij = 0. Let us now make a qualitative 
discussion on various contributions to rji and U. ViVj . 

The one-loop diagrams contributing to (or IT^o) can be divided into [DD], [DF] and 
[FF] types which involve two gauge vertices, one gauge and one Yukawa vertices, and two 
Yukawa vertices [TJ. The explicit formulae for Ud (and also IT n ) are presented in Appendix 
B. The [DD] type contains diagrams with Z/W, sneutrino/Higgs, charged slepton/ Higgs ex- 
changes. Among them, the Z/W exchanging diagrams give corrections of Oij 2 -^^-) ^ 10~ 4 
which are aligned with the tree level masses and thus can be safely ignored. The contribution 
from the charged slepton/Higgs exchange diagrams (involving also charged lepton/chargino 
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exchange) can be as large as ^tp apart from the possible alignment as can be read from 
mixing matrices ( |T2i|T3|) and (|A22| , [A23|) . The same is true for the sneutrino/Higgs exchang- 
ing diagrams (involving also neutrino/neutralino exchange). From our explicit calculation 
in the next section, we will see that the seemingly large correction of order j^tp for large 
tan/3 is compensated by a large suppression factor and thus 3> fy. The [DF] type in- 
cludes the diagrams with h T coupling involving charged slepton/Higgs exchange, A coupling 
involving stau/tau exchange, and A' coupling involving sbottom/bottom exchange. The 
[DF] diagrams with h T coupling can be negligible for small tan (3 as the correction is of the 
order < 1<T 5 and are sub-leading to the corresponding diagrams contributing to 

n n for large tan/?. In the [FF] type, there are diagrams with h T h T , h T X and h^X' couplings. 
Generically, these diagrams are sub-leading to the similar ones in U n . 

There are also [DD], [DF], and [FF] types in the IT n (or H uu ) diagrams. The [DD] type 
can be ignored as the Z exchanging diagrams give aligned corrections of the order of 1 % 
and the sneutrino exchanging ones are also negligible. In the [DF] type, there are diagrams 
with the coupling h T and A. The first ones give rise to the correction ~ which has 

tan/? enhancement compared to the corresponding diagrams in Up as mentioned before. In 
the [FF] type, there are diagrams involving h T h T , h T X, XX and A'A' couplings. The first two 
classes of diagrams may have significant contributions for large tan j3 as noted in Refs. [BT 



which will be partly confirmed in the next section. The last two classes of diagrams are 
the frequently considered ones as they typically give the main contributions to the one-loop 
correction. 

Our qualitative discussion on the one-loop corrected neutrino masses will be confirmed 
by numerical computation performed in the next section. For our explicit computation, 



we use gauge- mediated supersymmetry breaking models [17| where the universality arises 



more naturally than in supergravity models. A rigorous calculation in minimal supergravity 



models with bilinear terms has been performed in Ref. [12 



IV. ATMOSPHERIC AND SOLAR NEUTRINO OSCILLATIONS IN 
GAUGE-MEDIATED SUPERSYMMETRY BREAKING MODELS 

In this section, we discuss the properties of the one-loop improved neutrino mass matrix 
obtained in the previous section through explicit calculation. Then, our emphasis is put on 
addressing whether or how the atmospheric and solar neutrino masses and mixing angles can 
be realized in the context of R-parity violating models with supersymmetry broken by the 
gauge mediation mechanism. In a gauge mediation model, soft parameters of sfermion fields 
depend only on their gauge quantum numbers upon supersymmetry breaking, and thus the 



supersymmetric flavor problem is resolved in a natural way |T8[ . However, a difficulty in this 
type of models is that the origin of \x and B terms cannot be provided in its own context. 
This also applies to the R-parity violating bilinear terms /ij, B^. Assuming, however, B and 
Bi terms come from the same origin, we take B/Bi = Then, one can rotate away the 

R-parity violating bilinear terms into the Higgs bilinear term at the mediation scale and 
thus no sneutrino VEVs can be generated. But this property is spoiled by renormalization 
group evolution and thus sneutrinos get nonzero VEVs at the weak scale. We solve renor- 
malization group equation to determine these sneutrino VEVs for given input parameters 
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or A', A. In the R-parity preserving case, an extensive study of sparticle spectrum through 



renormalization group evolution has been done in Ref. ||19|| . In this work, we take a minimal 
gauge mediation model assuming one pair of messenger multiplets (5 + 5) |[7J and put the 
messenger scale equal to the messenger sector supersymmetry breaking scale A. In most 
cases, we will set A = 75 TeV. 

The simplest R-parity violating model is the "bilinear model" in which only bilinear 
R-parity violating terms are present at the fundamental scale, or at the mediation scale 
of supersymmetry breaking in our case. This model having only three R-parity violating 
parameters is more predictive and turns out to be very restrictive as we will show. It is also 
conceivable that there are both bilinear and trilinear couplings originated from a certain 
flavor symmetry which also explains quark and lepton mass hierarchies ||13|| . In this spirit, 



we assume that A« = Aj33 and A- = A^ 33 are much larger than the other couplings and 
give rise to experimentally relevant contributions. We call this the "trilinear model" which 
has five R-parity violating parameters (note that A3 = 0), as the bilinear terms /ij can be 
rotated away into the trilinear couplings. For our purpose, it is then convenient to use 
renormalization group equations in the basis where the bilinear terms fii are rotated away 
at all the scales which are recollected in Appendix C. 

Before investigating realistic solutions to neutrino problems, it is worth discussing how 
much one-loop corrections can change tree-level results. As an example, we take a bilinear 
model with = /i 2 = {13 = 10~ 4 /i at the mediation scale A = 75 TeV. By rotating away /ij, 
we have 

K = ^h h} Xi — —h T . (24) 

In Table I, we show the values of R-parity violating parameters and one-loop contribu- 
tions to sneutrino VEVs calculated at Q 2 = m^m^. In the table, q, <p~, S°, P°, ip~, 
and Z/W denote the squark, charged slepton/Higgs, CP-even and CP-odd sneutrino/Higgs, 
neutrino/neutralino, charged lepton/chargino and gauge boson contributions to E^, respec- 
tively. The main contribution to E^ comes from stops or sbottoms, and we do not show 
vanishing or negligible contributions from ip° or bottom quark, etc. For small tan/5, E^ 
gives about 18% correction to the tree- level quantity, B{t^ + rriL i H 1 [see Eq. @]. The dom- 
inant contribution to Y,^ comes from charged slepton/Higgs. Contrary to the case of E^, 
the squark contribution is negligible as there is an almost precise cancellation between up- 
type and down- type squark contributions. The correction E^ 2 ) to m 2 is typically small and 
5% in our case. Taking both one-loop contributions E^ 1 ' 2 ' into account, the one-loop im- 
proved sneutrino VEVs, £j, get increased by 13 % from the tree-level value for tan/3 = 3. 
This corresponds to a 27% increase of the tree-level neutrino mass as it is proportional to 
£ 2 . This is a significant change when one tries to explain the experimental data. More 
seriously, the one-loop correction could even upset the tree-level calculation. As can be seen 
from Table I with tan/? = 30, the one-loop improved value £j is about 7 times larger than 
the tree-level value £f. This phenomenon occurs due to a cancellation between the tree- level 
quantities. As one can see from Table I, even though each quantity rn 2 L . Hi or Bjtp is about 
10 times larger than the one-loop correction E^, there exists a cancellation of first two 
digits in the sum, Bit/3 + m^H^. Therefore, it is essential to include one-loop corrections in 
the calculation of neutrino masses and mixing arising from R-parity violation. 
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To show the tan (3 dependence of tree and loop masses, we present order of magnitudes 
of some mass matrix components in the last five rows of Table I. Let us first note that while 
rji get increased mildly with tan /3, the components of n n with h T \i couplings get increased 
tremendously with tan /3 and become comparable to the contributions from AA couplings, as 
was first noted in Ref. pi| . While the heaviest neutrino mass dominated by the tree- level 



value gets little changed, the second largest neutrino mass strongly depends on tan (3. 

Let us now turn to an important question whether a realistic neutrino mass matrix can 
be obtained in our scheme. We first discuss the bilinear model. As can be seen in Table 
I, a bilinear model with low tan/5 produces two heavy neutrino masses in the right range 
for the atmospheric and solar neutrino vacuum oscillation mass scale. But, it appears to be 
impossible to obtain bi-large mixing required in this case. In order to have bi-large mixing, 
one needs £/j 3 <C 1 which means £i <C £2,3 as the tree-level mass dominates [[J. This implies, 
in turn, /ii <C /i 2 ,3 which gives rise to the solar neutrino mixing angle ~ /ii//U 2 ,3, leading to a 
contradiction. On the other hand, one can easily find a small mixing MSW solution to solar 
neutrino problem |22| for a reasonably large tan (3. A solution is as follows. Let us take the 



input values with tan/3 = 24; 

^1 = 6x10-6, ^ = ^ = i.5xl(T 4 . (25) 

This gives rise to the neutrino mass matrix in the eV unit; 

/5.7xl(r 5 1.4xl0~ 3 1.3xl(T 3 \ /2.2 x KT 5 5.5 x KT 4 2.3 x 10" 4 ' 

- 1.4 x 1(T 3 3.6 x 1(T 2 3.2 x 1(T 2 + 5.5 x 1(T 4 1.4 x 1(T 2 5.8 x 1(T 3 

Vl.3xl(T 3 3.2xl0- 2 2.9xl0- 2 / \2.3 x 10~ 4 5.8 x 10~ 3 2.5 x 10~ 3 

/3.5xHT 5 8.7xl0~ 4 l.lxl0~ 3 \ 
= - 8.7 x KT 4 2.2 x KT 2 2.6 x 10~ 2 , (26) 
\l.lxHT 3 2.6 x 10-2 2.6 x 10-2/ 

where the first matrix on the left-hand side is the tree-level mass and the second one is the 
one-loop mass. This leads to the mass eigenvalues and oscillation amplitude; 

mt = 4.6 x 10~ 19 , m 2 = 2.3 x 10" 3 , m 3 = 5.1 x 10~ 2 

sin 2 26 sol = 0.0035, sin 2 26 atm = 0.99, sin 2 26 chooz = 0.0029 . (27) 

In our notation, sin 2 29 so \ is the z/ e disappearance oscillation amplitude corresponding to 
Anijj, and sin 2 26* atm (sin 2 26*^007) is the — > v T {y e disappearance) oscillation amplitudes 
corresponding to Am| 2 31 , which are measured in the solar neutrino 0, and the atmospheric 
neutrino [J| (CHOOZ 0) experiments, respectively. Let us note that such a solution can 
be obtained in a limited range of tan/3 = 22 — 26. If one reduces the mediation scale, a 
larger tan/3 is needed to achieve this kind of solution; e.g., tan/3 ~ 40 for A = 50 TeV. Our 
investigation shows that the small mixing MSW solution to the solar neutrino problem can 



be obtained for rather smaller tan/3 than claimed in Ref. [21|. The main reason for this 



discrepancy is that the possibility of a precise cancellation for large tan/3 discussed above 



was not considered legitimately in a qualitative analysis of Ref. [21]. In fact, for tan/3 > 30 



(with A = 75 TeV), the loop corrections get weird due to a large third generation effect 
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1) so that the input, /i2 
does not work at all. 



/13, motivated by large atmospheric neutrino mixing, 



In trilinear models, we can have more freedom and thus find more variety of solutions. 
A remarkable result is that the bi-maximal mixing solution for the atmospheric and solar 
neutrino vacuum oscillations [23| can be obtained in our framework. For tan/3 = 2, a set of 
the input values to provide bi-maximal mixing is 



X 1 = 2.3 x 1(T 6 , 
A; = -9.4 x KT 7 



Ao = 7.7 x 1(T 6 



An — Aq 



-1.3 x 10 



-5 



(28) 



A key feature to find such a solution is that one needs X[ smaller than A' 2 3 by about one 
order of magnitude to yield <C 1 and rather large Ai to increase the mass components 
in the v e direction. The resultant mass matrix is given by 



/1.3 x 10" 4 
1.8 x 1(T 3 

\ 1.7 x iir 3 



1.8xl0- 3 1.7 xlO" 3 
2.7 xlO" 2 2.5 x 10- 2 
2.5 x 10~ 2 2.3 x 10~ 2 
1.2 x 10~ 4 
1.8 x 10~ 3 




6.7 x 10~ 6 
2.9 x 10~ 5 
1.6 x 10~ 5 



2.9 x 10~ 5 
1.1 x 10~ 4 
5.0 x 10~ 5 



1.6 x 10~ 5 
5.0 x 10~ 5 
-2.2 x 10~ e 



x 10~ 3 1.7xl0~ 3 ' 



2.7 x 10~ 2 2.5 x 10~ 2 
1.7xl0~ 3 2.5 x 10~ 2 2.3 x 10~ 2 

leading to the following mass eigenvalues and oscillation amplitudes; 

m! = 2.4xl0- 6 , m 2 = 8.1xl0" 6 , m 3 = 5.0 x 10~ 2 
sin 2 26 sol = 0.97, sin 2 2fl atm = 0.99, sin 2 2d chooz = 0.0098 



(29) 



(30) 



In this case, we have Am^ 



0.6 x 10 _10 eV 2 . Compared to the other solutions, this is 
rather sensitive to the input values, but one does not need to rely on fine-tuning of the input 
parameters. 

Our last example is to accommodate the small mixing MSW solar neutrino solution even 
for low tan/3. It is the case where Aj is larger than A^. Specifically, the input values with 
tan {3 — 5; 



Ai 

a; 



1.6 x 10~ 6 
4.9 x 10~ 7 



pve rise to the mass matrix 

'4.0 x KT 5 
1.0 x 10~ 3 



1.0 x lO" 3 8.5 x 10" 4 
2.5 x 10" 2 2.1 x 10~ 2 
8.5 x 10- 4 2.1 x 10" 2 1.8 x 10~ 2 



A 2 = 4.0 x 10~ 5 
Ai = A', = 1.2 x 10" 



/ 4.4 x 10~ 6 
1.1 x 10~ 4 
V -7.5 x 10~ e 



(31) 



1.1 x 10~ 4 
2.7 x 10~ 3 



X. iu ±.0 A iu j y — I .0 A iu — J 

/3.6xl0- 5 8.9 x 10~ 4 8.5 x 10- 4 

= - 8.9 x 10~ 4 2.2 x 10~ 2 2.1 x 10~ 2 

1 ' x 10~ 4 2.1 x 10~ 2 1.8 x 10~ 2 



-7.5 x 10~ 6 ' 
-1.9 x 10~ 4 
1.9 x 10~ 4 -1.4 x 10~ 6 



(32) 



The corresponding mass eigenvalues and oscillation amplitudes are 



12 



mi = 1.5 x 1(T 20 , m 2 = 1.4 x 1(T 3 , m 3 = 4.2 x 10~ 2 

sin 2 26 sol = 0.0029, sin 2 26 atm = 0.99, sin 2 2# chooz = 0.0035 . (33) 

In this solution, the loop mass matrix is dominated by the diagrams with AA couplings and 
the other diagrams can be neglected. In all the realistic solutions, we find that the mass 
matrix elements coming from the diagrams with h T X and h^X' in IT^, and with gh T and h T h T 
in Ho,n are smaller than the others by more than one order of magnitude. Furthermore, 
the Ud diagrams can be omitted as they give negligible contributions to the resulting mass 
eigenvalues and mixing angles. Among the IT n diagrams, one should not ignore the diagrams 
with gX, h T X in the parameter region with large tan/3, or relatively large A (as in the first 
two solutions), as their contributions become comparable to those with AA or A'A' which 
usually give largest contributions. 

V. CONCLUSIONS 

We have performed a comprehensive analysis on the one-loop corrected neutrino masses 
and mixing. Our first concern was to determine sneutrino vacuum expectation values at 
the one-loop level by employing one-loop effective scalar potential. R-parity violation mixes 
R-parity even and odd particles and thus sfermion and fermion mass matrices are more 
complicated than in the R-parity conservation case. However, under the condition that 
R-parity violating couplings are small as implied by small neutrino masses, one can make 
an analytic diagonalization up to the quadratic order in R-parity violating parameters to 
compute the complete one-loop contributions to sneutrino minimization condition. The 
results are presented in Appendix A. This one-loop determination of sneutrino VEVs is 
essential to confront the R-parity violating models with experimental data. It has been 
observed that the one-loop correction can give rise to a significant change from the tree-level 
calculation, in particular, when there occurs cancellation between tree-level contributions. 
Taking gauge-mediated supersymmetry breaking models, we have shown that this kind of 
cancellation arises naturally for large tan j3 and the one-loop correction can change even the 
order of magnitude of the heaviest neutrino mass. 

We have then considered one-loop renormalization of neutrino mass matrix to determine 
mass eigenvalues and mixing angles for all three neutrinos. Concerning the neutrino sector 
only, it turns out that the use of the Q-independent effective mixing matrix is very useful 
to calculate one-loop corrected neutrino masses and mixing. In this procedure, we compute 
7x7 one-loop neutrino/neutralino mass matrix in the weak basis, and make a Q-independent 
seesaw diagonalization to obtain 3x3 effective neutrino mass matrix. From this effective 
mass matrix, one can simply calculate mass eigenvalues and mixing matrix working at zero 
momentum p 2 = 0. It is then argued that this procedure simplifies the calculation in the 
basis of fj,i — if one wants to calculate the leading contributions which are summarized in 
Appendix B. 

Following the above procedures, we have investigated the properties of neutrino mass 
matrix to see whether the atmospheric and solar neutrino oscillations can be obtained in the 
framework of R-parity violating theories with gauge mediated supersymmetry breaking. In 
the bilinear model, together with the atmospheric neutrino oscillation, one can only obtain 
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the small mixing MSW solution to the solar neutrino problem taking large tan (3 ^ 20. 
For the case of the trilinear model, one can have more solutions. With low tan/5, the bi- 
maximal mixing solution to the atmospheric and solar neutrino problems can be realized. 
Furthermore, varying 5 free parameters, it has been shown that the small mixing MSW 
solution can be obtained even for small tan j3. We would like to emphasize that all of these 
R-parity violating solutions to neutrino problems can be realized without any fine-tuning, 
apart from the smallness of R-parity violating parameters, /ij, A and A', which may be 
originated from a certain flavor symmetry explaining the quark and lepton mass hierarchies. 



APPENDIX A: ONE-LOOP CONTRIBUTIONS TO SNEUTRINO VEVS 

We collect here one-loop contributions to sneutrino minimization condition ([!]) by calcu- 
lating all the field-dependent particle masses. 

Stop contribution: The stop mass matrix including R-parity violating contribution is 
given by 

ka2 = fm 2 Q3 + m\ + |(4M^ - Mf)c 2/3 A t v 2 + h t {^v x + mm) \ ( ., 
1 \ A t v 2 + htifivi + fitUi) m 2 iR + m\ -\{M&- M 2 )c w ) 

where h t and A t are the top quark Yukawa coupling and trilinear soft parameter respectively, 
and M^c 2 p = g 2 w (v \ + u 2 — v f) , M z c 2 p = g 2 (v 2 + u 2 — v 2 ). Following the procedure described 
in Section II, we find the contributions from the stop mass eigenstates ti >2 as follows: 

6 I AM 2 Dt m t fii 




{ j) ~ 64tt 2 1 2 v 2 + 2m? - m? - m 2 6v 2 ' ^ I ^ 2 ' 1 ' 



3 64:% 2 } 2m 2 — m? — m? v 2 s 2 rI2 

6 flM 2 M 2 ,-M 2 Rt s.Ur, -hMl 




where M 2 t = M 2 U , M|, t = M 2 12 and M\ t = M 2 22 ignoring small R-parity violating 
contributions. 

Sbottom contribution: The sbottom mass matrix is 



U 2_ ( m Q s + < ~ U 2M w + M 2 z )c 2 p -A b v 2 -h bf i Vl + A' iUi \ 

J b ~ { ~A b v 2 - hifJLV! + A[ Ul m\ + m 2 + |(M^ - M 2 )c w ) [K6) 



The sbottom contributions are 




in 2 I 1 V\ 2m? — m 

'2',/ , 6 [ lMj M 2 Lb -M 2 Rb -AMl + M 2 

u 1 jJ 64vr 2 I 2 2m- 2 - m? - m- 2 6w 2 

^ 6w bi b 2 




Egfo) = -— { 2^— + - 2 T ^ )ml\-^-l\ (A4) 

1 J 647T 2 1 ui 2m? - m? - m? n ^ [ b ' \ Q 2 " 




where M| 6 , Mr&, and M|, 6 are defined as in the stop case. The contributions from the first 
two squark generations can be obtained by obvious substitutions. 
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Charged slepton/Higgs contribution: Let us denote the charged slepton and Higgs con- 
tributions to S's as 

,2 \ 



where stands for the mass eigenstates en, and H . Here Sj (<f>) = {drrA/ duj) /2v\ and 

Sj 2 \<p) = (dm^/duj)/2uj can be calculated from the diagonalization of the mass matrix ( 
H|6D. In the matrix the R-parity conserving mass components are 

Ml = ml + mf - \{2M 2 W - M 2 z )c w , 
M 2 Di = -A\v x - h\v 2l i 

M 2 Ri = ml iR + mf + {M 2 W - M 2 )c w (A6) 

apart from u 2 dependence in the gauge boson masses. The slepton mass diagonalization 
matrix elements c e j and s e « define the rotation angle as tan 6^ = s e i/c e i and useful quantities 
which are free from the ambiguity in defining the angle 9 e i are 

sin2g CT = 2 1Md \ , cos26 ei = M ? M f (A7) 
mi — m- mi — m~ 

A straightforward calculation leads to 
Sf\en, So) = {2<« 1 A i F 1 - (m% Hl + B^c'F, 

+(m 2 LiHl t p - B t ) ({M^ - m\.)Mlels w + M*^ - m^M^) 

X* + A*/?)^ + (m^ i )(m»(Ml i - m|_)/cj 



i i 
mi — m~ 

en e t 2 



+(ml iHl tf3 - (m 2 , - m 2 H -)Mlc 2 s 2p + M%^ 
~'" D ' [A^Fi + mlmtpFi 



m- — mi 

en ei2 



—(rntfn) ({M 2 Li + M 2 m - 2m\-){M 2 u - m 2 H - + 2M^^)tp - 2m^M 2 Di /c 2 
+krn 2 LiHl tp ~ Bt) (m^{M 2 Ll + M 2 m - 2m\-) - M 2 Dl {M 2 Ll - m 2 H - + 2M^)s 2p )] } 



2 

S?\H-) = - <g^_g! ( (A ^ - m^M^faw + M^^ + M^m> 

i ^A 1 * f^e../i\s2 _2 \ / 3 /jiyr2 _2 , oA/f2 „2 < 



(m^{M 2 Ll - ro|_)/cj - M^(M^ - m 2 H . + 2M 2 w c%)t^ (A8) 



where Fi = (mf^ — m 2 H -){m\. 2 — m 2 H S). 

Note that the contributions for the two light sleptons are reduced to the following simple 
form taking m\ 2 = 0: 
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Similarly, one finds 



rajgi + Bjtp m 2 L . Hl t p - Bj M^s w 



(M 



Li 



m 



(A9) 



Sf\eii 1 e i2 ) 



I Ml 



4 v 2 



± 



M 



Li 



mi — mi 



'3Mf 
4 w 2 



M 2 W 



+ 



1 A 



2 ufFi 



{(m 2 ^ 2 - Mi, - M^4) F 4 + mf /i 2 (M 2 , - m 2 H _)/c 



+ (M 2 t - m 2 H -)(M 2 Li - m\. + 2M^4) 2 S J - 2m^t^M 2 Dl (M 2 u - m 2 H - + 2M 2 w c 2 p) 



M 2 Ll 



M 2 Rl 



~ m li 



2M 2 , 



(m%-4 - Ml - M 2 W 4) Fi - mfn 2 (Mh - m 2 H -)/c 2 



+ (M 2 m - m 2 H - ) {M 2 Li - m 2 H . + 2M^c 2 ) 2 S 2 



2A lVl - M 2 Di ) Fi - M 2 Dl (mff, 2 /4 + {M 2 Ll - m 2 H - + 2M 2 r c 2 ) 2 s 2 ) 
+ mt^Ml + M 2 m - 2m 2 H .)(M 2 Ll - m 2 H . + 2M^c 2 )] } 



S?\H-) = ^ { (m 2 w {1 + 2sl) - \m 2 z c w + (M 2 Lt - m^)t 2 ) Fi 



(M 2 4 - m^)(Ml - m^ + 2M^)% + mf ^{M 2 Li - m^)/c 2 



+2m^t f3 M 2 Dl (Ml - m\„ + 2M^cJ)} 
CP-even sneutrino/Higgs contribution: As in the above case, we write 

s£ 2) (0) 



2 SP>(0)m 2 (ln^-l 



64vr 2 



(A10) 



(All) 



where <p runs for Re(z/), /i° and H° and S"} (0) are calculated from the matrix ( |I~4"D as 
follows: 

{m 2 LiHi (m 2 A0 c 2 p - m?) + Bitp{m A0 c 2 p + m?.)' 



Sj 1} (*<) 



S^\h°,H ) = ^^-^{(mi^^oCa^-TOy + Bitp(m%c 20 + ml.] 

Z V Pi 

Tc 2a (m 2 LtHi (m 2 A0 cl + Mf - mf.) + ^(m^os 2 + Mf - m?.; 

T S -f ((™W/» + S0(m^ + Mf )c w + 2(m 2 iiHl ^ - £?>)<)} 
Mf(m? -m 2 c^) > 



S?\h°,H ) 



lMf 
2 

1M| 
2 w 2 



1 + ^11 + 



2 

mi. 



Fi 

"F c 2a 



1 + 



Mz(m 2 A0 - ml )c 2 p 



TS2aS 2 



Mfm? 



(A12) 
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where Fi = m~. — 2m|.(m^ + m 2 H0 ) + m\ m 2 H0) and m~. = m|. + ^M z c 2 /3- 
Recall that 



_ //t A o - ivi z _ ,,u a q -r iKzg . 

C2a — — C2/3 2 2~> S2a ~~ ~ s 2/3 2 T~ l A1<3 J 

CP-odd sneutrino/Higgs contribution: Again, one can calculate from the matrix (p 



E & 2, W = ^P'(0K(lng-l), (AM) 
where runs for Im(i>), A . Here, s!- ^ 2 \(j)) are given by 

gj 1) (i>/) = -^ < V^ ' ^ (1) ^ = 

Gauge boson contribution: 

Bottom quark contribution: 

sg ) (0) = O (A17) 

N eutrino /neutralino contribution: To diagonalize the 7x7 neutrino/neutralino mass ma- 
trix, we use the seesaw formula to obtain the rotation matrix: 

where the seesaw diagonalization matrix separating out the neutrino and neutralino 
parts is given by Gtv = MpM^ 1 , and U u is the neutrino mixing matrix which is determined 
from the full neutrino mass matrices including the tree and one-loop contributions. The 
matrix N is the usual neutralino diagonalization matrix whose analytic expression can be 
found, e.g., in the paper by Barger et al. in Ref. |]L4}|. We present a simpler form as follows: 



N 2l 


1 m v o — Mi 


Nu 


t w m x o — M 2 


N 3l 


1 (m x o - Mi)(m x oc /3 - ixsp) 


N u 


s w M z {m x o - (is 2 p) 


N 3l 


1 (m x o - Mi)(m x osp - ixcp) 


N u 


s w M z (m x o - fis 2f3 ) 



(A19) 
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and the Nij satisfies the unitarity relation, + N 2i + + = 1. Recall that one has 
to keep R-parity violating contributions to calculate S's which we do not show explicitly in, 
e.g., flM9|). 

The sneutrino VEVs contributions are 



^(1,2) 



8 



64vr 2 



Q' 2 



(A20) 



where ip runs for four neutralinos Xii an d S^'^ty) are given by 

M 2 

S\ X) {^) = m—fui^ \-m\tp + m^M^tp + n) - /iMJ /D(^) 



^ 2) (V0 



Ml 



2 ?7ty 



-mj + m^M^ + m^ 2 - ^ilfJ /Dty) , 



(A21) 



where 



Dty) = 5mJ - Aml(M l + M 2 ) + Sm^M^ - /i 2 - M 2 Z ) 
+2m^[/i 2 (M 1 + M 2 ) + M 2 z {2^s 2 p + Mj)] + det(M 



Note that S!- 1 ' = in the basis with /ij = 0. 

Charged lepton/chargino contribution: Diagonalization matrices U R ,U L bring respec- 
tively the weak eigenstates (e°, W + , H 2 ) and (e«, W/ - , ifi ) to the mass eigenstates (ef , xf 2 ) 
and are given by 



U 



R 



( 1 



—r 



Cy S y 

/ 1 l\Ci + li,Si -l\Si+l\Ci 



V-4 



Sz 



where r^i and Zi >2 are given by 

fii V2mlfj 1 M w (fj 1 - M 2 tp)c[3 



■ Ui_ 



/x (M 2 /x + M 2 ,^) 2 
,Ui _ jjj V2jjM w C/3 



n 



^1 



(A22) 



//' (M 2 /i + M 2 ,^) 2 



/i M 2 /i + M^s 2/3 /i 



/i M 2/ u + M^s 2/3 

and the useful combination of the usual chargino mixing angles are given by 
M 2 m- + urn 

Xi 



+ 



(A23) 



SlS r 



X 2 



2 2 

m _ — m - 

Xi x 2 
urn - + M 2 m 

A1 



X-2 



2 2 

m _ — m 



ClS r 



SiC r 



\/2M w (spm x - + cpm x -) 



2 2 
m _ — m 



Xi x 2 

y/2M w (cpm - + spin 



X-2 



x 2 



2 2 
777. _ — 777 



(A24) 



Xi 



x 2 
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Neglecting the contributions from the two light charged leptons by taking m 



Sp^m 2 I In- 



- 1 



64vr 2 1 vr/ V\ Q 2 , 
where ip runs for r and the charginos X12 an d S^ 1,2 \ip) are given by 



0, we get 
(A25) 



Sf >(V0 = /i^ [2M^(// - M 2 ^)(mJ - m 2 ) + fem^ (2M^ + (mj - M 2 2 )/c 2 )] /D^) 



— Ajm T — - 



2 

x3 



(A26) 



Sf } (V>) = 2^ [mjM^ - m 2 M^Gu 2 + 2il^ S 2 + m 2 ) + m 2 M^(/i 2 + 2M^sj 

+5 i3 m 2 T {m£/2cg + m^(M^(3 - ^) - Mf/cg) - 2M 2 v (M 2f j,t (3 - 3M^s| - /i 2 )}] / D{ip) 
where 



-Df'0) = 5mf,, — 4mi(m 2 - + m 2 _ + m 2 ) + 3m 2 ,,[m 2 (m 2 - + m 2 _) + m 2 -m 2 -l — 2m 2 m 2 _m 2 - . 

vry V ^ v Xi X 2 T ' ^ L rV Xi X 2 7 Xi X 2 J r Xi x 2 

APPENDIX B: ONE-LOOP CORRECTIONS TO NEUTRINO MASSES 

Here we summarize the calculation of one-loop self-energy diagrams in the weak basis. 
A key ingredient in this procedure is to express the propagators inside loops in terms of 
the mass basis propagators which can be done by the use of the diagonalization matrices 



15, 


18, 


A18 


A22) 



propagators to the mass basis ones as follows: 



(Bl) 



where a and (3 represent weak basis fields and ip mass eigenstates. For neutrino/neutralino, 
charged lepton/chargino, sneutrino/Higgs (scalar or pseudoscalar) , and charged slep- 
ton/Higgs propagators, one finds 



1 P>°2 



7J N U N 



v 4>° _ +h,S,P ? ,S,P 



Vt + =U L - I M R +I + 
p 1 pZ Pi ^ Plr 

<p<t>° _ hjs,p 7J s,p 



<p4 



<Pi Pi 



u 



c 



-U C +th+ 

tpj<t>+ 



vB2) 



where p's and ?/>'s represent the neutrino/neutralino or charged lepton/chargino fields in the 
weak and mass basis, respectively; (p° is a complex sneutrino or Higgs field and 0° is a scalar 
or pseudoscalar mass eigenstates; <p~ and 0~ represent the charged slepton/Higgs fields in 
the weak and mass basis, respectively. 

We first present Il n which gives most significant contributions to the one-loop neutrino 
self-energy. 
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^ E ^"- e c^ i+ ^- S o(p 2 ,^-,mJ-) 



*is E <- r c(^-<e 7 J^-5o(p 2 ,mJ_,mJ_) 



16vr 2 

ij>-,4>- 

Y^l E K^4^ m ^- fi o(p 2 ,™ 2 -,mJ-) 

y v „ 

~^ E m ^o(p 2 , m 2 , m? ) + (i <- j) (B3) 



167T 2 Y ^ ' ^ 

The contributions to Ud are 



i/j ,4, 

9 hi tw \ - r_-py>~ -p</>~ _ -p^~ t>*~ 
16vr 2 V2 , . 6l< e ^ H i + # i" e * 

-n^^ E p*;[2P&j,-p&ii^-Bo(P 2 ,m|., ro 5-) 

V if) ,(p 

-^^ep^U+pUw^x.^) (B4) 

g 2 1 

-^4 E [-^>C^]<^o(p 2 ,<,0 
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9_K i 
16tt 2 ^2 



T.[-HU^ m2 Mp 2 ^ m l m l) ( B5 ) 



K 



1 

tp° ,(jfi 

_9ht_ r_T*T I?* - -, +T>^~ -p^~ J m 2 Enfr? 2 m 2 m 2 ) 

V>-,0- 

E h^^ 1+ -^ e ^J<5o(/,m 2 _,<) 



16tt 2 



l>!X Ty^^Y^B^,mi,mi) (B6) 



16vr 2 Y bxLbiR 

g 2 1 X - r _,/,0 _. U _, ,'l 



^2 



167T 2 ,^ L ' W-W+- e iL H- 

-jt* E [^I^^-K-BbCp'.rnJ-.mJ-) (B7) 
Note that is equal to 1 for ^~ = and vanishes for all others, and that TV's. = TCbSb, 

1 i L R 

V^jp = c 2 , s 2 , and T^g* = s 2 , c 2 where c& and s& define the sbottom diagonalization angle; 
tan 6*b = Sb/cb as in Eq. (&7). The function £> is defined by 



i? (p 2 ,m 2 ,m 2 ) = ^L^ln4-ln^ + l with p 2 - 0. (B8) 



APPENDIX C: RGE IN THE BASIS WITH NO BILINEAR TERMS 

Renormalization group equations for the lepton number violating parameters in the basis 
where the bilinear terms /ijLjifi in the superpotential are rotated away have been calculated 
in the second paper of Ref . pO] . Here we collect the results with a few minor corrections. 



= KiSisK + h\ + 3h 2 - lg< 2 - 3gl - fgj) (CI) 

16tt 2 ^ = UZhl - 3g' 2 - 3gl) (C2) 

(7 A' 7 1fi 

16tt 2 — i = + h\ + 9/> 2 - V - 3fe 2 - A 2 ) + MZhK) (C3) 

at 9 3 

7 16 

+ 2A<($ i3 /i T A r + Mt + 3/^6 + -^Mi + 3(? 2 M 2 + — ^Ma) 
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16tt 2 — - = M5h 2 T - 3g' 2 - 3g 2 2 ) + A'^hK) + A,(6/i r A r + 6g' 2 M 1 + 6g 2 2 M 2 ) (C4) 

16tt 2 — |^ = m\ iHl (5 l3 h 2 T + h 2 T + 3h 2 b ) - 6A[A b - 2A t A T (C5) 
- 6\%(m 2 L . + m 2 Q3 + m 2 D3 ) - 2\ih T (m 2 L . + m| 3 + m| 3 ) 
16tt 2 ^ = ^(fe/i 2 + 3h\ - g' 2 - 3g 2 2 ) - fi(6A% + 2A l h T ) (C6) 

Note that these equations are valid when the R-parity violating couplings are small enough 
so that the next leading corrections are negligible. 



22 



REFERENCES 



[1] L. Hall and M. Suzuki, Nucl. Phys. B231 (1984) 419. 

[2] A. S. Joshipura and M. Nowakowski, Phys. Rev. D51, 2421 (1995); M. Nowakowski and 

A. Pilaftsis, Nucl. Phys. B461, 19 (1996); F. M. Borzumati, Y. Grossman, E. Nardi 

and Y. Nir, Phys. Lett. B384 (1996) 123; B. de Carlos and P. L. White, Phys. Rev. 

D54 (1996) 3427; A. Yu. Smirnov and F. Vissani, Nucl. Phys. B460 (1996) 37; R. 

Hempfling, Nucl. Phys B478 (1996) 3; H. P. Nilles and N. Polonsky, Nucl. Phys. B484 

(1997) 33; E. Nardi, Phys. Rev. D55 (1997) 5772. 
[3] E. J. Chun, S. K. Kang, C. W. Kim and U. W. Lee, Nucl. Phys. B544, 89 (1999), 

|hep-ph/ 9807321 ; A. S. Joshipura and S. K. Vempati, |hep-ph/ 98082321 [hep-ph/ 9903431 ; 

J. Ferrandis, |hep-ph/9810371| ; M. Bisset, O. Kong, C. Macesanu and L. Orr, [hep 



ph/9811498| ; S. Rakshit, G. Bhattacharyya, Phys. Rev. D59, 091701 (1999); R. Ad- 



hikari, G. Omanovic, Phys. Rev. D59, 073003 (1999); D. E. Kaplan and A. Nelson, 
hep-ph/990l251 ; A. Abada and M. Losada, |hep-ph/ 9908352] ; J. L. Chkareuli, et al, 
|hep-ph/ 9908451 ; J. C. Romao, et al, |hep-ph/9907499| ; O. Haug, et al, |hep-ph/99093T8| . 

[4] Super-Kamiokande Collaboration, Y. Fukuda, et.ai, Phys. Rev. Lett. 81, 1562 (1998). 

[5] B. T. Cleveland et al, Astrophys. J. 496, 505 (1998); Kamiokande Collaboration, K. 
S. Hirata et al, Phys. Rev. Lett. 77, 1683 (1996); GALLEX Collaboration, W. Hampel 
et al, Phys. Lett. B447, 127 (1999); SAGE Collaboration, J. N. Abdurashitov et al, 
astro-ph/9907113| Super-Kamiokande Collaboration, Y. Fukuda, et al, Phys. Rev. Lett. 
82, 2430 (1999). 

[6] LSND Collaboration, C. Athanassopoulos et al, Phys. Rev. Lett. 81, 1774 (1998). 
[7] CHOOZ Collaboration, M. Apollonio et al, Phys. Lett. B420, 397 (1998); [Kep 
ex/99070371 



[8] B. Mukhopadhyaya, S. Roy and F. Vissani, Phys. Lett. B443, 191 (1998). 

[9] E. J. Chun and J. S. Lee, Phys. Rev. D60, 075006 (1999), |hep-ph/981120T| ; S. Y. Choi, 
E. J. Chun, S. K. Kang and J. S. Lee, Phys. Rev. D60, 075002 (1999), [hep-ph/ 9903461 . 
[10] S. Kiyoura, M. M. Nojiri, D. M. Pierce and Y. Yamada, Phys. Rev. D58, 075002 (1998). 
[11] R. Hempfling, in Ref. @. 
[12] J. C. Romao, et al, in Ref. [|]. 

[13] V. Ben-Hamo and Y. Nir, Phys. Lett. B339, 77 (1994); T. Banks, Y. Grossman, E. 
Nardi and Y. Nir, Phys. Rev. D52, 5319 (1995); E. J. Chun and A. Lukas, Phys. Lett. 
B387, 99 (1996); K. Choi, et al, Phys. Lett. B394, 89 (1997); Phys. Rev. D60, 031301 
(1999); G. Eyal and Y. Nir, JHEP 9906, 025 (1999), |hep-ph/99044T3; ; G. K. Leontaris 
and J. Rizos, |hep-ph/9909206| ; Q. Shafi and Z. Tavartkiladze, |hep-ph/ 9909238 . 



[14] R. Arnowitt and P. Nath, Phys. Rev. D46, 3981 (1992); V. Barger, M. Berger and P. 

Ohmann, Phys. Rev. D49, 4098 (1994). 
[15] M. Nowakowski and A. Pilaftsis, in Ref. f|; 

[16] D. Pierce and A. Papadopoulos, Nucl. Phys. B430, 278 (1994); Phys. Rev. D50, 565 
(1995). 

[17] M. Dine and A. Nelson, Phys. Rev. D48, 1277 (1993); M. Dine, A. Nelson and Y. 

Shirman, Phys. Rev. D51, 1362 (1995); M. Dine, A. Nelson, Y. Nir and Y. Shirman, 

Phys. Rev. D53, 1658 (1996). 
[18] For a review, see G. F. Giudice and R. Rattazzi, [hep-ph/9801271 . 



23 



[19] S. Dimopoulos, S. Thomas and J. D. Wells, Nucl. Phys. B488, 39 (1997). 
[20] B. de Carlos and P. L. White, in Ref. @; E. J. Chun, el al, in Ref. @. 
[21] K. Choi, K. Hwang and E. J. Chun, Phys. Rev. D60, 031301 (1999); D. E. Kaplan and 
A. Nelson, in Ref. @. 

[22] For a recent analysis, see J. N. Bahcall, P. I. Krastev and A. Yu. Smirnov, Phys. Rev. 
D58, 096016 (1998). 

[23] V. Barger et al, Phys. Lett. B437, 107 (1998); H. Georgi and S. L. Glashow, [hepj 
| ph/ 98082931 ; A. J. Baltz et al, Phys. Rev. Lett. 81, 5730 (1998); S. Davidson and S. 
F. King, Phys. Lett. B445, 191 (1998); C. Giunti, Phys. Rev. D59, 077301 (1999); 
Y. Nomura and T. Yanagida, Phys. Rev. D59, 017303 (1999); R. N. Mphapatra and 
S. Nussinov, [hep-ph/98083011 ; M. Tanimoto, |hep-ph/98075~T7] ; H. Fritzsch and Z. Xing, 
Phys. Lett. B440, 313 (1998); S. K. Kang and C. S. Kim, Phys. Rev. D59, 091302 
(1999); C. S. Kim and J. D. Kim, |hep-ph/9908435| ; C. Jarlskog et al, Phys. Lett. B449, 
240 (1999); R. Barieri et al, |hep-ph/990T24T[ ; Q.Shafi and A. Tavartkiladze, Phys. Lett. 
B451, 129 (1999); A. H. Guth et al, JHEP 9908, 018 (1999); J. Ellis and S. Lola, Phys. 
Lett. B458, 310 (1999); M. Jezabek and Y. Sumino, Phys. Lett. B457, 139 (1999); J. 
A. Casas et al, Nucl. Phys. B556, 3 (1999); Y. -L. Wu, |hep-ph/990522l D. Falcone, 
hep-ph/ 990531^ ; G. AltareUi et al, |hep-ph/ 9905536) ; A. Ghosal, |hep-ph/990547q ; C. H. 
Albright and S. M. Barr, |hep-ph/9906297| . 



24 



TABLES 



TABLE I. R-parity violating parameters m 2 L Hl ,Bi, and and one-loop corrections in the 
GeV unit. Each collum corresponds to the lepton generation, i = 1, 2, 3. are one-loop improved 
sneutrino VEVs and ^ are tree values in the unit of v\. The last five rows show typical sizes of 
some components of tree and loop mass neutrino matrices and their eigenvalues in the eV unit. 





tan (3 = 3 


tan = 30 


Bi tan p 


4.00 x IO" 2 
-3.95 x IO" 3 


4.00 x IO" 2 
-3.95 x IO" 3 


3.93 x IO" 2 
-3.88 x IO" 3 


3.67 
-3.60 


3.67 
-3.60 


3.59 
-3.52 




4.93 x IO" 3 


4.93 x IO" 3 


4.93 x IO" 3 


5.02 x 10" 1 


5.02 x 10" 1 


5.00 x IO" 1 


</>- 


1.39 x 10~ 3 


1.39 x 10~ 3 


1.45 x 10~ 3 


-1.73 x 10~ 2 


-1.73 x 10~ 2 


-3.89 x 10~ 2 


s° 


-5.64 x 10~ 4 


-5.64 x 10~ 4 


-5.02 x 10~ 4 


-5.03 x 10~ 2 


-5.03 x 10~ 2 


-4.61 x 10~ 2 


dO 


7.21 x 10~ 4 


7.21 x 10~ 4 


7.55 x 10~ 4 


1.48 x 10~ 3 


1.48 x 10~ 3 


1.39 x IO -3 


_,fll 


6.47 x IO" 3 


6.47 x IO" 3 


6.63 x IO" 3 


4.36 x 10" 1 


4.36 x 10" 1 


4.17 x IO" 1 


o 

m i>i 


57280 


57280 


62860 


56440 


56440 


61340 


Q 


-4.84 


-4.84 


-4.84 


-6.22 


-6.22 


-6.22 


<T 


2507 


2507 


2867 


2382 


2382 


2178 


S° 


-356.3 


-356.3 


-360.2 


-329.8 


-329.8 


-331.9 


pO 


739.6 


739.6 


935.6 


759.3 


759.3 


928.9 




-712.6 


-712.6 


-712.6 


-486.3 


-486.3 


-486.3 




782.8 


782.8 


787.3 


745.0 


745.0 


1048 


z/w 


-258.8 


-258.8 


-258.8 


-258.8 


-258.8 


-258.8 




2696 


2696 


3253 


2805 


2805 


3072 




-7.09 x IO" 7 


-7.09 x IO" 7 


-6.35 x IO" 7 


-8.61 x IO" 6 


-8.61 x IO" 6 


-7.50 x IO" 6 




-6.30 x 10" 7 


-6.30 x IO" 7 


-5.63 x IO" 7 


-1.32 x IO" 6 


-1.32 x IO" 6 


-1.08 x IO" 6 


tree mass 




2 x 10~ 3 






3 x IO" 3 




n^, : 77 




10~ 5 






io~ 4 




n n : /iA 




10~ 8 






io- 2 




n n : AA 




IO" 5 






10~ 2 




eigenvalues 


o, 


1 x 10~ 5 , 


7 x 10~ 3 


o, 


7 x IO" 3 , 


2 x 10~ 2 
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